Critical aging of Ising ferromagnets relaxing from an ordered state 
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Abstract 



We investigate the nonequilibrium behavior of the d-dimensional Ising model 
Cj ■ with purely dissipative dynamics during its critical relaxation from a mag- 

I . netized initial configuration. The universal scaling forms of the two-time re- 

^ , sponse and correlation functions of the magnetization are derived within the 

O \ field-theoretical approach and the associated scaling functions are computed 

up to first order in the e-expansion (e = A — d). Aging behavior is clearly 
displayed and the associated universal fluctuation-dissipation ratio tends to 
^ . X^ = |[1 — (^ — ^)e + 0{e^)\ for long times. These results are confirmed 

^^ \ by Monte Carlo simulations of the two-dimensional Ising model with Glauber 

tH- \ dynamics, from which we find ^^^c ~ 0.73(1). The crossover to the case of re- 

'nI" ■ laxation from a disordered state is discussed and the crossover function for the 

\^ ' fluctuation-dissipation ratio is computed within the Gaussian approximation. 

o' 

g : I. INTRODUCTION 

'^ ' The dynamics of statistical systems close to critical points has been the subject of in- 

5 ■ tensive theoretical and experimental investigations in the last three decades, during which 

O ■ mainly equilibrium properties have been studied in detail. As in the case of static properties, 

the presence of a nearby critical point greatly facilitate the study of dynamical behavior. In 
K^ \ fact, upon approaching a second-order phase transition some of the relevant features of the 

;h ' dynamics (at length and time scales much larger than the microscopic ones) are determined 

only by the fluctuations of a suitable order parameter and they become to some extent in- 
dependent of the actual microscopic details of the system {universality) . In turn, this allows 
a quantitative characterization of the dynamical behavior within the various universality 
classes grouping together the microscopically different systems which display the same crit- 
ical behavior. This can be usually done to an extent that is in general well out of reach for 
the general case of non-critical dynamics. 



Critical dynamics provides also a simple instance of slow collective evolution. Statistical 
systems with slow dynamics have recently attracted considerable theoretical and experimen- 
tal interest, in view of the rich scenario of phenomena they display: Dramatic slowing down 
of relaxation processes, hysteresis, memory effects, etc. After a perturbation a system with 
slow dynamics does not generically achieve equilibrium even for long times and its dynam- 
ics is neither invariant under time translations nor under time reversal, as it should be in 
thermal equilibrium. During this never-ending relaxation aging occurs: Two-time quantities 
such as response and correlation functions depend on the two times s and t > s not via t — s 
only and their decays as functions of t are slower for larger s. At variance with one-time 
quantities (e.g., the order parameter) - converging to asymptotic values in the long-time 
limit - two-time quantities clearly bear the signature of aging. 

Aging was known to occur in disordered and complex systems (see, e.g., Ref. [1]) and 
only in the last ten years attention has been focused on simpler systems such as critical 
ones, whose universal features can be rather easily investigated by using different methods 
and which might provide insight into more general cases [2] (see also Ref. [3] for a pedagog- 
ical introduction). Indeed, consider a system with a critical point at temperature Tc, order 
parameter (^(x, t), and prepare it in some initial configuration (which might correspond to 
an equilibrium state at a given temperature Tq). At time t = bring the system in contact 
with a thermal bath of temperature Ti,{^ Tq). The ensuing relaxation process is expected 
to be characterized by some equilibration time teq(Tfe) such that for t ^ tcq(^b) equilibrium 
is attained and the dynamics is stationary and invariant under time reversal, whereas for 
< t ^ teq(Tf,), the evolution depends on the specific initial condition. Upon approach- 
ing the critical point Tf, = T^ the equilibration time diverges and therefore equilibrium is 
never achieved. To monitor the time evolution we consider the average order parameter 
M{t) = ((/9(x, t)) (translational invariance in space is assumed), the time- dependent corre- 
lation function of the order parameter C^it, s) = ((/^(x, t)ip{0, s)), where (. . .) stands for the 
mean over the stochastic dynamics, and the linear response (susceptibility) Rx{t,s) to an 
external field. R^it, s) is defined by -Rx(^5 s) = 6{ip{'x,t))/6h{s), where h is a. small external 
field, conjugate to (/^(x = 0, s) (e.g., if y? is a magnetic order parameter, h is the magnetic 
field), applied at time s > at the point x = 0. Note that causality implies R^it, s > t) = 
and that Cx{t,s) = Cx{s,t) in the bulk. According to the general picture of the relaxation 
process, one expects that for t > s ^ teq(Tft), Cx(t, s) = C^(t — s) and Rx{t, s) = R'^it — s) 
where C^ and R^ are the corresponding equihbrium quantities, related by the fluctuation- 
dissipation theorem (FDT) 

1 dC^Hr) 
K'ir > 0) = -JTdy- ' ^^^ 

where the temperature is measured in unit of the Boltzmann's constant. The FDT suggests 
the definition of the so-called fluctuation- dissipation ratio (FDR) [4,5]: 

where we assume t > s. According to the previous picture for t > s ^ tcq(Th) the FDT 
yields X^{t, s) = 1. This is not generically true in the aging regime. The asymptotic value 
of the FDR 



X°^ = lim limXx=o(t,s) (3) 



+00 1 



is a very useful quantity in the description of systems with slow dynamics, since X°° = 1 
whenever the aging evolution is interrupted and the system crosses over to equilibrium 
dynamics, i.e., teq(Tf,) < 00. Conversely X°° 7^ 1 is a signal of an asymptotic non-equilibrium 
dynamics. Moreover X°° can be used to define an effective non-equilibrium temperature 
Teff = T/X°°, which might have some features of the temperature of an equilibrium system, 
e.g., controlling the direction of heat flows and acting as a criterion for thermalization [6]. 
Within the field-theoretical approach to critical dynamics it is more convenient to focus on 
the behavior of observables in momentum space. Accordingly hereafter we mainly consider 
the momentum- dependent response Rq(t, s) and correlation Cq(t, s) functions, defined as 
the Fourier transform of Ryc{t, s) and Cy^{t, s), respectively. In momentum space it is natural 
to introduce a quantity that, just like Xx(t, s), "gauges" the distance from equilibrium 
evolution [7]: 

Note that A:'q(t, s) is not the Fourier transform of X^(t, s). The long-time limit 

^„^o= lim limA'q=o(t,s) (5) 

has been used to define an effective temperature, in analogy with X°°. Interestingly, it has 
been argued that X°° = X^q [7] (see also Ref. [8]). However, for quenches to the critical 
point X^Q (and therefore X°°) depends on the specific choice of the quantity which R^ and 
Cq refer to [10] and therefore it would be difficult to assign a sound physical meaning to T^s 
defined from X°° (see also [9]). 

We shall consider here the universal aspects of one of the simplest dynamic universality 
class: Pure relaxation (Model A in the notion of Ref. [11]) within the Ising (static) univer- 
sality class. As we shall discuss later on, the lattice Ising model with Glauber dynamics 
belongs to it. In spite of its simplicity the non-equilibrium behavior of this model has been 
investigated in some detail only recently. In particular, earlier works focussed primarily on 
general scaling properties of the non-linear relaxation [12,13] close to T^ and on the equation 
of motion describing the dynamics of the order parameter in the presence of time-dependent 
external fields [14,15]. Two-time quantities such as correlation and response functions were 
not discussed explicitly beyond equilibrium [14], although the non-equilibrium behavior of 
Rq=o{t,s) was implicitly encoded in the equation of motion computed in Refs. [14,15]. A 
more careful analysis of the relaxation process from a given initial (non-equilibrium) state 
revealed the presence of an additional stage of relaxation with universal features [16,17] 
(characterized by the initial-slip exponent 6, cf. Sec. Ill) which was previously overlooked. 
The general analysis of Refs. [16,17] provided the scaling functions for generic multi-time 
quantities and in particular for the order parameter M(t) and for the two-time quantities 
R^{t,s) and C^{t,s). On the other hand there was no explicit analysis of the resulting 
aging behavior at the critical point, which has been later carried out in Ref. [5,7,18-20,10] 
by computing the two-time response and correlation functions for the case of a relaxation 
from an initially disordered state with vanishing order parameter (corresponding to a quench 



from high temperature to the critical point). In particular it has been shown that X°° = 
for quenches below T^, whereas X°° is a universal amplitude ratio for quenches to T^, as 
pointed out by Godreche and Luck [18]. Accordingly it is possible to take advantage of 
this universality by computing X°° within the field-theoretical approach to critical dynam- 
ics. For the Ising universality class with purely relaxational dynamics this analysis leads 
to X°^{d = 3) = 0.429(6), X°°{d = 2) = 0.30(5) in good agreement with the Monte Carlo 
results X^cid = 3) ^ 0.40 and X^^id = 2) = 0.340(5) (see Ref. [2] for details). These 
investigations have been also extended to different dynamic universality classes both in the 
bulk [2,21-23,8,24,25] and at surfaces [26,2,27]. 

In this paper we address the complementary problem of the relaxation occurring at the 
critical point from an initial state with non-vanishing value of the order parameter. Taking 
advantage of previous results (mainly reported in Refs. [17,15]) we provide predictions for 
the universal scaling forms of two-time response and correlation functions, discussing the 
resulting universal aspects of aging behavior within the field-theoretical approach to critical 
phenomena. [Within this approach some equal-time correlation functions were previously 
investigated: Cx(t, t) in Ref. [17] whereas C(^=o(t,t) in Ref. [28], in the context of a detailed 
study of the relaxation in finite volume - see also Ref. [29].] We also compare the provided 
predictions with the results of Monte Carlo simulations of the two-dimensional Ising model 
with Glauber dynamics. The same problem has been recently addressed for other universal- 
ity classes [30,9,31]. The comparison of our findings with these recent results is presented 
in the conclusive section. 

The paper is organized as follows. In Sec. II we introduce the model and its basic 
features. In Sec. Ill we derive the general scaling forms for the correlation and response 
functions of the order parameter after a quench from an ordered state. In Sec. IV we 
solve the model in the Gaussian approximation. Such an approximation allows us to fully 
describe the crossover from ordered to disordered quench. Then in Sec. V we present our 
first-order e-expansion computation for the response and the correlation functions and from 
them we calculate the FDR. Sec. VI is devoted to an extensive Monte Carlo simulation of 
the two-dimensional Ising model with Glauber dynamics. Finally in Sec. VII we summarize 
our results and compare them with the existing literature. 

II. THE MODEL 

The simplest non-trivial model in which aging occurs is probably the Ising model in d 
dimensions evolving with a purely dissipative dynamics after a quench to the critical point. 
Its Hamiltonian on an hypercubic lattice is given by 

n = -Y.s.s^, (6) 

(xy> 

where S^ = ±1 is a spin located at the lattice site x. The sum runs over all pairs (xy) 
of nearest-neighbor lattice sites. A purely dissipative dynamics for this model proceeds by 
elementary moves that amount to randomly flipping the spin S'x. The transition rates can 
be arbitrarily chosen provided that the detailed-balance condition is satisfied. For analytical 



studies the most suited are the Glauber ones [32], which allow exact solutions in the one- 
dimensional case [33-37,24]. 

Despite its simplicity, the dynamics of this model is not exactly solvable in physical 
dimensions d = 2,3. In order to investigate analytically the dynamical behavior we take 
advantage of the universality considering the time evolution of a scalar field (/)(x, t) (ob- 
tained in principle by coarse-graining the lattice spin variables Sx) with a purely dissipative 
dynamics (Model A in the notion of Ref. [11]). This is described by the stochastic Langevin 
equation 



()V9(x,t) 



where fl is the kinetic coefficient, ^(x, t) a zero-mean stochastic Gaussian noise with 

(e(x,t)e(x',t')) = 2f]5(x-x')5(t-t'), 



(7) 



(8) 



and Ti-lip] is the static Hamiltonian. Near the critical point, 7i[(/9] may be assumed of the 
Landau-Ginzburg form 



^M 



d'^x 



-(Vv?)' + -roV^' + -W 



(9) 



where tq is a parameter that has to be tuned to a critical value tqc in order to approach the 
critical temperature T = T^ {tqc = 0, within the analytical approach discussed below), and 
go > is the bare coupling constant of the theory. This coarse-grained continuum dynamics 
is in the same universality class as the lattice Ising model with spin-flip dynamics [11]. 

Correlation and response functions of a field satisfying the Langevin equation (7) can be 
obtained by means of the field-theoretical action [38,39] 



S[ip, ^] 



r 
dt / d'^X 



dip 



ipVtip 



(10) 



where (^(x, t) is an auxiliary field, conjugate to the external field h in such a way that 
7i[(/9, h] = Hlip] — J d'^xhip. In terms of S, the average over the stochastic dynamics induced 
by the noise ^(x, t) of a quantity 0[p\ depending on the order parameter field p is given 

by [38,39] 



{0[p)]) = [[dpdip] 0[p] e-^l'^-'^] 



(11) 



The effect of the external field h on this average (denoted, for /i 7^ 0, by (. . .)h) is accounted 
for by using 7i[v?; h] in Eq. (10). As a consequence, the linear response to the field /i of a 
generic observable O is given by 



5{0)h 



5h{:x., s) 



fi(^(x,s)0) 



(12) 



h=0 



and hence "^(x, t) is termed response field. In particular 



i?.-.(M)-'^^^"''^^^ 



(5/i(x', s) 



= fi(^(x',s)^(x,t)). (13) 

h=0 



The effect of a macroscopic initial condition ^poi'^) = vi^^t = 0) may be accounted for 
by averaging over tlie initial configuration with a weight e"^''^'^"] where [16] 

H,[ipo] = ld'x^[M^)-Mo]^ (14) 

which specifies an initial state with Gaussian short-range correlations of the order parame- 
ter fluctuations, proportional to Tq-^ and spatially constant initial averaged order parameter 
Mq = {(po{x.)) (the generalization to space-dependent Mq is straightforward). In the ex- 
perimental protocols we have in mind the system has initially a non-vanishing value of the 
order parameter (magnetization). This can be obtained by preparing the system either (a) 
in an equilibrium low-temperature state in the absence of external fields or (b) by applying 
an external field at generic temperature. In both cases the resulting state can be described 
by Eq. (9) with external field h, leading, far enough from the critical point, to a Gaussian 
distribution such as Eq. (14) with tq ~ QoMq where Mq = {—Qro/goY^'^ in case (a) (ro < 0) 
and Mq = (6/1/(^0)^ in case (b). However, within the renormalization-group (RG) approach 
to the problem it has been shown [16] that Tq"^ is an irrelevant variable in the sense that 
Tq-^ affects only the correction to the leading long-time scaling behavior we are interested 
in. In view of that we fix it to the value Tq"^ = from the very beginning of the calculation. 
Accordingly, cases (a) and (b) previously mentioned differ only in the corrections to the 
leading scaling behavior. 

In order to account for a non-vanishing mean value of the order parameter ((/^(x, t)) = 
M{t) (we assume that M{t) stays homogeneous in space after the quench) during the time 
evolution, it is convenient [15] to write the action (10) in terms of fluctuations around M(t), 



^(x,t) = v?(x,t)-M(t), V^(x,t) = ^(x,t), (15) 

so that ('?/'(x, t)) = {ip has been introduced to make the notation uniform). The problem 
we shall consider is therefore the dynamics of fluctuations of the fleld '?/'(x, t) in the time- 
dependent "background fleld" provided by M{t). For later convenience we also introduce 
the rescaled magnetization m{t): 

2 M^ 

m =9o—, (16) 

so that a perturbative expansion in go leads to a flnite value for rriQ = m(t = 0) in the case 
(a) mentioned above. The resulting action in terms of i/JjI/j may be written as 

/•OO f 

S = J dtj d^'xiCo + £1 + £2) , (17) 

with 



Cq = -^pn^p + ij(dt + n [v^ + ro + m^{t)] ij) 



9o /j.\ 7 / 2 , ^6*0 7 , 3 



2 



C2 = i^{dt + n 



ro + :^rri^it) 




—m{t) =^hcsit). 
90 



(18) 
(19) 

(20) 



We split up the action S so that Co is the Gaussian part, £i contains the interaction 
vertices and £2 gives the couphng to the effective magnetic field hcs(t) acting onip{q= 0, t) 
and due to a nonzero M{t) oc m{t). Note that in Co the effect of a non- vanishing mean 
value of the order parameter m{t) is equivalent to a time-dependent temperature shift: 
(m(t),ro) ^-» {0,ro + m'^{t)). Therefore, even if the system is asymptotically (for long times) 
at its critical point tq = tq^ = 0, for finite times it is effectively in the disordered phase. 

Following standard methods [38,39] the response and correlation functions may be ob- 
tained by a perturbative expansion of the functional weight e~^^^'-^''^^~^^°^'-^°^^ in terms of the 
coupling constant go- The propagators (Gaussian two-point functions of the fields ip and ip 
in momentum space) are [16] 



(V^(q, t)V^(q', t'))o = (27r)^5(q + q')<(t, f) , 
(^(q, t)^(q', t'))o = (27r)^5(q + q')C°(t, t') , 



with (g = |q|) 



i?°(t, t') = 0{t - t') exp l-n (g2 + ro)(t - t') + f dt"m^{t 

fOO 

C°(t, t') = 2Q dt"Rl{t, t")Rl{t', t"). 



(21) 
(22) 



(23) 
(24) 



In the following we will assume the Ito prescription (see, e.g., Refs. [17,38]) to deal with 
the ambiguities arising in formal manipulations of stochastic equations. Consequently, all 
the diagrams with loops of response propagators have to be omitted. This ensures that 
causality holds in the perturbative expansion [16,17,39]. From the technical point of view, 
the breaking of time-translation invariance does not allow the factorization of connected 
correlation functions in terms of one-particle irreducible ones as usually happens when time- 
translation invariance holds. As a consequence, as in the case of surface critical phenomena 
[40], the whole calculation has to be done in terms of connected functions only [16]. 

The perturbative expansion can be as usual organized in terms of Feynman diagrams 
with propagators given by Eqs. (23) (represented as a directed line with the arrow pointing 
to t, the larger of the two times t, t') and (24) (represented as an unoriented line) and 
vertices given by Ci. In addition to the standard time-independent quartic vertex —^go a 
time-dependent cubic one —Qy/2gQm{t) , due to a non- zero magnetization m{t), has to be 
accounted for in the expansion. As explained below, the contribution of hcs{t) cancels for a 
suitable choice of m{t). The evolution equation for the magnetization m{t) can be obtained 
by solving the equation of motion {6S/6ilj{:K,t)) = 0. Taking into account that causality 
implies (■?/'(x, t)) = (corresponding to — 4A = 0) and that m{t) has been defined so that 
(■?/'(x, t)) = 0, one finds (in zero external field) [15] 




1 
ro + -m^(t) 



) m(t) + n^mit) (^2) + ^ (V^^^) = , (25) 



which has to be supplemented by the proper initial condition, e.g., by assigning the value 
of tuq = m{t = 0) oc Mq. Note that the sum of the last two terms in this equation 
equals —Ti, where 7^ = ^^^ is the one-point vertex due to the interaction vertices in £i, 
whereas the first term is —7^, where 7^ = -^-U is the contribution of the effective field 
hesit) [see Eq. (20)] coming from £2- Accordingly, Eq. (25) states that 7^ + 7^ = 0, i.e., 
^^t(5) = -<\i) + ^^-n = where ^^^ is the total one-point vertex with one external ^-\eg 
due to £1 + £2- In view of this simplification we will omit, in the following computation, all 
the diagrams having the one-point vertex as a subdiagram. To the lowest order [hereafter 
(dq) = d'^q/{2TiY and Vt = I] 



(^^(x, t)) = J {dq)C'^{t, t) + 0{go) , (^^(x, t)) = O(^) (26) 

and therefore up to one loop and at the critical point Tq = 0, Eq. (25) becomes 



= dtm{t) + -nx\t) + -f- I {dq)C% t) + 0{g',) . (27) 



III. SCALING FORMS 

The non-equilibrium evolution of a critical system in the absence of an order-parameter 
background is a well-understood topic (see, e.g., Ref. [2]). In particular, it is well-known 
that zero-momentum response and correlation functions satisfy the scaling forms [16,2] 

i?q=o(t, s) =AR{t- sTit/sYj^Ris/t) , (28) 

Cq=o(t, s) = Ac sit - snt/sfj^cis/t) , (29) 

where a = [2 — 1] — z)/z, z is the dynamical critical exponent, r] the anomalous dimension 
of the field, and ^ is a genuine non-equilibrium exponent [16]. Ar and Ac are non-universal 
amplitudes which are fixed by the condition J^r,c(0) = 1. With this normalization J-r^c 
are universal. From these scaling forms the universality of X°° follows as an amplitude 
ratio. Although in what follows we focus mainly on the case q = 0, the generalization of 
the presented scaling forms to non- vanishing q amounts to the introduction of an additional 
scaling variable y = A^Qq^lt — s). Aq is a dimensional non-universal constant which can be 
fixed according to some specified condition. 

These behaviors are clearly changed if a critical system evolves starting from a state 
with Mo 7^ 0. It is known from general scaling and RG arguments that the magnetization 
satisfies the scaling form [16] 

m{t) = Arnm^e+'j^MiBmm^t'') (30) 

where k = 6 + a + f3/{vz) = 9 + j35/{vz) and standard notation for critical exponents 
has been used (note that /t is generically expected to be positive). The universal scaling 

8 



function J-'m{v) has an analytic expansion around v = 0, and, due to the obvious symmetry 
(mo,m(t)) h-i> (— mo,— m(t)), Tm{.—v) = J-'m{v). Accordingly the non-universal amplitudes 
Am and Bm can be determined, e.g., by imposing J-'m{0) = 1 and J-'m{0) = —1. From 
Eq. (30) one sees that the effect of a non- vanishing initial magnetization is the introduction 
of an additional macroscopic time scale Tm into the problem, i.e., Tm = {Bmrno)'^^'^ and of an 
additional associated scaling variable u = t/Tm in the scaling forms. In the long-time limit 
M 3> 1 one has to recover the well-known behavior of the critical relaxation m{t) ~ t-Pi^^ 
[13], and therefore J-'m{v ^ 1) ~ v"^. In that limit 

m{t) = ^t-/5/-^^,^((B„mo)^/'^t) (31) 

where the universal scahng function J^m,oo{u) = M.q + ^A^l u^'^ + 0{u^'^). 

Having understood the presence of a new scaling variable u = (Bm'moy^'^t = t/r^ 
(alternatively v = Bmmot'^ = m"), it is trivial to generalize the scaling forms for the response 
Eq. (28) and connected correlation functions Eq. (29) to the case of a non- vanishing mQ-. 

i?q=o(t, s) = Anit- snt/sfFnis/t, B^mon , (32) 

Cq=o(t, s) = Ac s{t - snt/sfFci-s/t, Bmmon , (33) 

where no new non-universal amplitudes have been introduced and obviously Fr{x^ 0) = 
J-'r{x) and Fc{x, 0) = J-'c{x). The resulting functions Fr and Fq are then universal. 

According to Eqs. (32) and (33) a non-vanishing mean value of the initial magnetization 
mo 7^ affects the scaling properties of the response and correlation function as soon as 
BrntTiot'^ ~ 1 (i.e., t ~ Tm) and in particular this happens in the long-time limit we are 
interested in, characterized hy t ^ s ^ Tm- This formally corresponds to the case mo — * oo, 
as opposed to the case previously considered, mo = 0. In this limit one expects the scaling 
forms (32) and (33) to turn into: 

i?q=o(t, s) = AR{t- sy{t/sf^R{s/t) , (34) 

Cq=o(t, s) = Ac s(t - snt/sffc{s/t) , (35) 

which resemble those for mo = 0. As before Ar^c are non-universal constants which are 
fixed by requiring ^r,c{0) = Ij where ^Rfi{x) are universal functions related to the large-u 
behavior of Fr^c{x,v). The "new" exponents 6 and 6 are clearly different from 6. In fact, 
for the related problem of the response and correlations of fluctuating modes transverse to 
the direction of the decaying magnetization in 0{N) models, it has been argued [31] that 
6 = 6= —f3/{vz). On the other hand, nothing is known, in general, about the longitudinal 
fluctuations that are the only degrees of freedom of the Ising model. In particular it is not 
obvious, a priori, whether 6 and 6 should be expected to be the same, if they are novel 
exponents - as ^ - or just combinations of known ones. The rest of this section is devoted 
to showing that Eqs. (34) and (35) hold for Model A with 






-6' = -^. (36) 

fZ 



In order to prove Eq. (36) we generalize the RG treatment of Refs. [16,17], confirming 
also Eqs. (32) and (33). The general scaling properties of the connected correlation function 
Gn,h,no = (M"'[<^]"'[<^o]"'°)c [where <^o(x) = "^(x, t = 0)] at the RG fixed-point and in the 
presence of a no n- vanishing initial magnetization (we restrict attention to the critical point, 
with Tq^ = 0) are given by [17] 

G'„,^,^„({x, t}; mo) = /^("'"'"«)G'„,^,,,,({/x, IH}; morn , (37) 

where 6{n, n, no) = n{d-2 + ri)/2 + h{d + 2 + fi)/2 + no{d + 2 + fi + ?7o)/2, {fj -rj)/2 = z -2, 
and fjo = —2z9. In order to determine the behavior of the correlation and response functions 
for s/t <^ 1 we perform a short-distance expansion in the time variable. Indeed, for s — > 0, 
one expects 

(^(x, s ^ 0) ~ (5(s,mo)(^o(x), (38) 

and therefore, inserting this relation into a generic correlation function Gn^n^ho-, one gets from 
Eq. (37) 

g(s,mo) = s-'Q(smy"), (39) 

where Q{x) is finite for x = [16], whereas the behavior for a; — ;► cxd has to be determined. 
According to the previous equation, 

VR^{t,s^ 0) = (v9(q,t)^(-q,s ^ 0)) ~ s-' Q{sml'''){^{ci,t)^o{-^)) , (40) 

where the volume V of the system has been introduced as a regularization of (27r)°'5(q = 0) 
following from translational invariance. Now we recall that (pq{(\ = 0) is the field conjugate 
to Mq oc rriQ (see Ref. [16]), thus 

(^(q = 0,«)^„(q = 0)) = ^M^^, (41) 

where, by definition, (v?(q = 0,t)) = VM{t) oc Vm{t), which scales according to Eqs. (30) 
and (31). Accordingly, the leading behavior in Eq. (41) is 

Smjt) j r+^ for t < r„, = (B^rrio)-^/- , 

6mo ^ \ t-i-'3/(-)mo"'-'/" for t > r^ , ^ ' 

Consider, first, the case s,t <^ r^ (corresponding to rriQ — > 0, i.e., r^ -^ oo). From the 
previous equations we recover the known result 

i?q=o(t, s ^ 0) ~ t" Q) (mo ^ 0) . (43) 

In the opposite hmit r^ ^ s,t (corresponding to mo -^ oo, i.e., r^ -^ 0) one gets Rq^=o(t, s — > 
0) ~ s^^ Q^suiq )t^^^f^^^''^^mQ ~ . Given that mo""*^ is irrelevant by dimensional analysis, 
we have that for mo -^ oo the response function must not depend on mo. Therefore Q{x — >■ 

oo) ~ x'^'^^ and 

10 



i?q=o(t, s ^ 0) ~ t« (^-j (mo ^ oo, i.e., mo > ^^^s^'/") • (44) 

According to this result we find for mo — > oo tlie scaling form Eq. (34) where 

«^-fl + o + A)._/!^. (45) 

V uz vz 



as anticipated. Note that by definition, we have 



F^{x, t; ^ oo) = ^x'-'^TKix) , (46) 

and so the ratio TZr = Ar/ Ar is universal, because it is the first term of the expansion of a 
universal function. 

Let us consider now the case of the correlation function, i.e., the short-distance expansion 
for the field (/^(x, s -^ 0). For mo = it is of the form (/^(x, s — *> 0) ~ Qi(s)(9sV5(x, s)|s=o, given 
that the Dirichlet boundary conditions implies v^(x, 0) = when inserted into correlation 
functions. Then one can prove diagrammatically that 0o(x) = (9s(/9(x, s)|s=o = 2(^o(x) when 
inserted into correlation functions. By means of this identity and of Eq. (37) one determines 
the scaling properties of Qi{s), leading to the expected scaling for the two-time correlation 
function (see Eq. (29)). In the present case, however, one has to take into account the 
non- vanishing initial value of the magnetization. Accordingly we expect an expansion of the 
form 

V9(x, s ^ 0) ~ Qo(s, mo)! + Qi{s, mo)ipo{x.) , (47) 

where 1 is the identity operator and the scaling form of the coefficients Qi can be determined 
by inserting Eq. (47) into Eq. (37). In particular this leads to Qo{s, mo) = s~^/*^'^^-'Qo(s^o )• 
As a consequence we have for the connected correlation function (the space dependence is 
understood) 

C{t, s^0) = im^s ^ 0)) = (v9(t)v9(s ^ 0)) - (v^(t))(v9(s ^ 0)) 

~ Qi(s,mo) Mt)^o) - {v{t)){^o)] ■ (48) 

Note that also in the case of non-vanishing initial magnetization the connected correlation 
function C^{t,t') [see Eq. (24)] satisfies the Dirichlet boundary condition [i.e., Cg{t,t' = 
0) = Cq(t = 0,t') = 0]. Together with causality this implies, as in the case mo = 0, that 
00 (x) = 2(^0 (x) when inserted into connected correlation functions, as can be verified by 
going through the diagrammatic proof presented in Ref. [16]. Therefore, by inserting (48) 
into correlation functions one finds 

gi(s,mo) = si-^Qi(smy'), (49) 

where Qi{x) is finite for x = 0. As before, the behavior for a; — > oo has to be determined. 
Taking into account that 
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VC^=o{t, s^0) = (v9(q = 0, t)v9(q = 0, s ^ 0)) ~ s^-'Silsm^'') {^{q = 0, t)^o(q = 0)) , 

(50) 

the reasoning proceeds as before (the only difference with Eq. (40) is the exponent of the 
prefactor s^~^), leading to the conclusion that Qi{x -^ oo) ~ x'^^^. Accordingly in the limit 
mo -^ oo one finds 

Cq=o(t, s) = Ac sit - snt/sffc{s/t) , (51) 

where 

Fcix, ^ ^ oo) = ^x'-'^cix) , (52) 

with universal ratio TZc = Ac/Aq- Interestingly enough the exponent 6 in Eq. (51) is the 
same as the one in Eq. (34), a result that could not have been inferred solely on the basis 
of scaling arguments applied to Eqs. (32) and (33). 

Comparing Eqs. (32) and (33) for mo = and Eqs. (34) and (35) (referring to the case 
mo -^ oo) one sees a fundamental difference between the quench from a disordered state and 
from the ordered one: In the latter case no novel exponent characterizes the aging properties 
and the non-equilibrium behavior is completely described in terms of equilibrium exponents. 

As a consequence of the scaling forms Eqs. (34) and (35) the limiting FDR can be written 
as 

which turns out to be a universal amplitude ratio related to X°°{m,o = 0). In Ref. [30] the 
question of the universality of the FDR for ferromagnetic models has been addressed and in 
particular X°°{m,Q ^ 0) has been computed for some effectively Gaussian models (see also 
Sec. IV). In particular, the fact that X°^(mo 7^ 0) 7^ X°^(mo = 0) has been interpreted as 
a signal of the existence of different dynamic universality classes for the "critical coarsening 
process" depending on mo being zero or not. Although it can be considered a semantic 
distinction, our point of view here is slightly different: X°°(mo 7^ 0) is as universal as 
X°°(mo = 0) in the sense that both do not depend on the microscopic details of the systems 
(i.e., on the specific realization of the universality class) but only on general properties such 
as symmetries, dimensionality, and conservation laws. Instead, the difference between these 
two quantities is due to the fact that the dynamic properties of a system belonging to a 
given universality class depend on the additional scaling variable u = t/Tm in a universal 
way. In turn, for long times, u can be either zero (r^ ~ ruo -^ 00, i.e., mo = 0) or 
infinite (rm finite, i.e., mo 7^ 0), making the specific value of rrio^ irrelevant (as expected 
from RG arguments). We will confirm this important feature explicitly within the Gaussian 
approximation in the following section. 

Note that the FDR could have been defined by using the non- connected correlation 
function Cq'^(t, s) instead of the connected one Cq(t, s), since in equilibrium both of them 
satisfy FDT. In this case we have 
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C^'it, s) = C^{t, s) + VM{t)M{s) , (54) 

where V is the volume of the system, leading to (t > s) 

C^Ut, s » r„) = Ac sit - sr{t/sf^c{s/t) + VCljr^'^'''\s-^l^-^-^'> , (55) 

with Cm = C^/ goY^'^ AmM-o / Bm [see Eq. (31)]. In a finite but large enough volume V 
(i.e., V ^ ^{tY ~ t'^^^) and in the aging limit t ^ s ^ Tm = (BmniQ)^^^'^, the second term 
dominates since j3/{vz) < {35 /{vz) — a. As a consequence the FDR takes the trivial value 
X-(mo^O)=0. 

IV. GAUSSIAN APPROXIMATION 

A first (oversimplified in physical dimensions d = 2, 3) description of the dynamics of the 
model just introduced is given by the Gaussian approximation of the action S [see Eq. (10)] in 
which one keeps only the linear and quadratic terms in the fields (■?/', ■?/'), neglecting additional 
anharmonic terms: S'g ['?/', "0] = /o°°dt/d'^x(£o + ^2)- This amounts to setting g^ = 0, but 
only after having performed the rescaling given by Eq. (16). In spite of its simplicity, the 
Gaussian approximation gives exact results above the upper critical dimension, equal to 4 
in Ising systems. Consequently the result of this section have to be equivalent to those 
obtained for the fully-connected lattice {d = 00) in Ref. [30]. 

For go = and at the critical point ro = 0, the equation of motion (27) simplifies to 

dtmcit) = --ml{t) , (56) 



3 



whose solution is 



-1 



This expression agrees with the scaling behavior (30) and (31), given that, within the Gaus- 
sian model, 6 = a = whereas t = |. The non-universal amplitudes appearing in Eq. (30) 

are Am = 1 and Bm = \/2/3 and therefore Tm = {Bmmo)^^^^ = 3/ {2ml), leading to 
mait) = mo(l + t/rm)~^/^. For t ^ Tm, Tncit) ~ (2t/3)^^/^ and therefore a nonzero rriQ^ 
affects only the corrections to the leading scaling behavior. In this sense ttiq^ is irrelevant 
for large times. 

Using Eqs. (23) and (57) the response function is (t > s) 

Rl{t,s) = e-^'i^-^)-!>t'mUt') = f-l±^Y\-e(t^s) ^ (58) 

\ t + Tfji / 

and the connected correlation function 

C% s)=2r dt'i?° (t, t')Rl{s, t') = 2 '7/r ^13/2 /' ^t'{t' + Tmfe'^''' . (59) 

Jo ^ ^ [{t + Tm){S + Tm)\-^'^ Jo 
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Note that R'^(t,s) < R'^{t,s)\r^=oo, i-e., the system evolving from a configuration with 
mo 7^ is stiffer than in the case ttiq = and also correlations are reduced, in agreement 
with the fact that the system is slightly displaced from its critical point. For q = one finds 

RU(i.s)=('-±^Y\ (60) 

^^=°^'' '^ = 2 lis + rSit + r J]3/2 • (61) 

Comparing this results with r^ = 0, with the scaling forms (34) and (35) we can identify 
z = 2, a = 0, 6 = —3/2 [see Eq. (45)] in agreement with standard mean-field exponents 
{6 = 3, u = P = 1/2, and t] = 0) and determine Aji = 1, Ac = |, ^r{x) = 1, and f'c{x) = 1. 
Taking into account the Gaussian expressions for the response and correlation function for 
rriQ = (see, e.g., Ref. [2]), one finds Aji = 1 and Ac = 2, and therefore TZji = 1 whereas 

Let us consider in more detail the fiuctuation-dissipation ratio. Using Eq. (23) one easily 
finds (9ti?°(t, t') = 5{t - t') - [q2 + ro + m|(t)]i?° (t, t'), and therefore, from Eq. (24), 

POO 

dsCl{t,s) = 2J^ dt'Rl{t,t')dMl{s,t') 

= 2i?° (t, s) - [q2 + ro + m^(s)]Cq(t, s) . (62) 

Accordingly, at the critical point (ro = 0) 

which, as one can see from Eqs. (58) and (59), actually depends on the scaling variables 
u = s/Tjn = 2sml/3 and y = |q|^s = q^s whereas it is independent of t, a typical property 
of the Gaussian approximation. The function X{u = 0,y), corresponding to the case of 
a quench from an high-temperature state, has already been investigated in Ref. [7]: For 
y = one has A?(0, 0) = |, whereas it monotonically increases towards the asymptotic value 
X[0,y -^ oo) = 1 as y increases. This suggests that the violation of the FDT for long 
times s -^ oo comes only from the zero mode q = 0, whereas all the other modes properly 
equilibrate. Let us consider the case mo -^ oo, i.e., u -^ oo. From Eqs. (58) and (59) one 
finds 

8y^ 

Xiu = oo, y) = ^-- , (64) 

^ '^' 9 - 12y + 6y2 + 8|/4 - 3e-22/(3 + 2y)' ^ ' 

which monotonically increases from the value | at y = to 1 for y ^ oo [see Fig. 1(a)]. 
As in the previous case, the fiuctuation-dissipation theorem is satisfied for long times by 
fiuctuations with q 7^ 0, whereas it is violated by the zero mode q = (i.e., by the global 
order parameter of the system), for which the FDR takes the value |. As anticipated, 
this behavior coincides with that one found from the exact solution of the dynamics of the 
d-dimensional Ising model with Glauber dynamics, in the limit of large d [cp. Eq. (47) of 
Ref. [30] with w ^ y to Eq. (64)]. 
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FIG. 1. Scaling function of the Gaussian FDR Xq{t, s) = X{u = s/Tm,y = q^s) for (a) niQ -^ oo 
as a function of y (see also Fig. 1 of Ref . [30] ) and (b) q = and finite ttiq as a function of u. 



Let us now focus on the FDR for the zero mode q 
u = s/Tm, given by [see Eqs. (60) and (61)] 



(i.e., y = 0), as a function of 



X{u,y 



°'^l 



-I 



l + u 



,-4 



(65) 



This quantity increases monotonically from the value ^ 



0, to the value | for u 



for u 

cxD. This implies that the FDR for the global order parameter (the total magnetization in 
magnetic systems) in the long-time limit approaches | whenever the mean value of the initial 
magnetization rriQ is non zero (and the value is indeed independent of the actual mo 7^ 0), as 
is the case for quenches from the ordered state, whereas the FDR is equal to | for mo = 0, 
i.e., for quenches from a disordered initial state. The plot of Xq=o{s,t) as a function of 
u = sJTm is reported in Fig. 1(b). 

In the case of a quench from the disordered phase it was instructive to look at the FDR 
in the real space x which has the same asymptotic value as A:'q(t, s) but a slower approach to 
it [2] . The expressions for the space-dependent Gaussian response and correlation functions 
can be worked out by Fourier transforming the corresponding equations (23) and (24). We 
do not report the result for the general case but we focus instead on the expressions for 
X = 0, i.e., on the autoresponse and autocorrelation functions {t> s): 



R\{t, s) 
C°(t,s) 



idq)Rl{t,s) 
{dq)C'(t,s) 










dt'(t + s-2t')"2(t' + r„) 



(66) 
(67) 



where Ka = (47r) 2 , and the expression (57) has been used. Introducing the scaling variables 
X = s/t and u = s/Tm, and the function Ia{x,u) = J^ df [l + a;(l —2v))]~^{l + uv)^, one finds 

dsC'Ait,s) , il-x)i 



Kit.s) 



:i+u) 



1 Id 

1 +-U 2 



x,u) + dxld.^{x,u 



(68) 
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Note that la{0,u) = [(1 + u)^ — l]/{4u) and therefore the previous expression becomes, for 
x = 0, 



(Js^ AV'^ • 



x=0 



A comparison with Eq. (65) shows that 

dsC\(t,s) 



u 



(69) 



R\{t,s) 



dsClit,. 



x=0 



i?o (t, s) 



(70) 



q=0 



i.e., in the hmit t ^ oo the FDR computed from the observables at x = (local) is the same 
as that one computed from the observables at q = (global), for every value of -u = s/Tm- 
This agrees with the heuristic argument of Ref. [7], which can be extended also to the case of 
non- vanishing initial magnetization. This has been explicitly checked in the case of infinite 
dimensionality in Ref. [30] . 

Consider now the case x = s/t -^ 1, i.e., the regime characterized by a short time 
difference (5t = t — s^s-^ oo. One expects on general grounds [2,41] that the system is in 
quasi-equilibrium with R\{t, s) = dsC^it, s), as one verifies explicitly: Indeed the correlation 
function can be written in terms of the scaling variables x and u as 



Note that for a > 1, Ia{x,u) 

C° (t, s 



u 



xld{x,u) 



(71) 



[2(a-l)]-i(l 



u 



[1 — x) "+-"^[1 + 0(1 — x)], and therefore 



-1 



(t_s)-2+i[l + 0(l-s/t)] 



(72) 



which is the equilibrium form of the two-point correlation function at the critical point 

C^/''\t - s) [we recall that C°("i)(t - s) = e^i'^'^^^/q^]. From Eq. (66) one easily finds 

R^it, s) = R} {t — s)[l + 0(1 — s/t)]. Note that these conclusions, leading to -R^(t, s) = 
(9sO° (t, s), are independent of the actual value of -u = s/rm and therefore of niQ. 

In the following we shall consider the thermoremanent magnetization defined as p{t, s) = 
Jq dt' RA{t,t'). According to Eq. (66) it can be written in the scaling form 



p°(t,s) 



i^.t-i+i 



1 + ^ 

X 



xJd{x,u) 



(73) 



where Ja{x,u) = Jq dv (1 —vx) "(1 + t>-u)2, with Ja{0,u) = |[(1 + -u)2 — l]/u and Ja{x,u) = 
(a — 1)~^(1 + u)2(l — x)^"+^[l + 0(1 — x)]. Therefore within the Gaussian approximation 
one finds: 



P'it,s) _1 . -Jji^^^) 

C%{t,s) 2^ ' n{x,u) 



(74) 



which gives, for x = 
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FIG. 2. Scaling function of p{t,s)/CAit,s) as a function of the scaling variables x = s/t for 
fixed u = s/Tm, within the Gaussian approximation for various values of the dimensionality d. 
(The Gaussian approximation becomes exact for d > 4.) From top to bottom u = cx), u = 1, and 
u = 0. Some of the qualitative features displayed by the Gaussian approximation are also found in 
the numerical results presented in Sec. VI (cp. the plot for d = 2 to Fig. 7(a)). 



p'{t,s) 



ciit, 



x=0 



4, ,3 1+M2 

-(1 + m)2^ {- 



(75) 



Note that, although for generic values of u this expression differs from i?^(t, s)/(9sC^(t, s) 
given in Eq. (69), they have the same asymptotic values for u = and u -^ oo. On the 
basis of scaling forms it is easy to realize that this equality holds beyond the Gaussian 
model and indeed in Sec. VI we shall compute the limiting FDR via the ratio between the 
thermoremanent magnetization p{t,s) and the correlation function CA{t,s). In Fig. 2 the 
ratio p°(t, s)/C^{t, s) is shown for fixed values of -u = s/Tm (from top to bottom: u = oo, 
u = 1, -u = 0) as a function of a; = s/t and of the dimensionality d. In spite of the 
fact that the Gaussian approximation provides exact results only for rf > 4 it displays 
correctly some qualitative features also for d < 4, although reliable quantitative predictions 
can be obtained only after having accounted for fluctuations. In particular Fig. 2 shows 
that p^{t, s)/C^{t, s) is an increasing function of x for fixed u, with an increasing concavity 
upon decreasing the dimensionality d. In particular the function is well approximated by a 
straight line followed by a sudden increase towards the value 1, for x — *> 1. This behavior 
is more pronounced upon decreasing d. The limiting value for x — > 1 is expected to be 1 
also beyond the Gaussian approximation, because it is a signature of the quasi-equilibrium 
regime, whereas the asymptotic value for x ^ depends not only on u (as correctly displayed 
by p°(t, s)/C^{t, s)) but also on the dimensionality d (we shall verify in the next section for 
the case u = oo). 



17 



V. ONE-LOOP FLUCTUATION-DISSIPATION RATIO 

In this section we present the one- loop (i.e., 0(e) in the e-expansion where e = A — d) 
perturbative computation of the correlation and response functions for a quench from a state 
with initial magnetization mo — ^ oo to the critical point. As we have argued from general 
scaling arguments and explicitly shown within the Gaussian approximation, a finite ttiq only 
gives corrections to the leading long-time behavior. The first step in this derivation is the one- 
loop equation of motion obtained in Refs. [14,15] which we rederive in the next subsection 
for the sake of completeness. Then we calculate the connected correlation and response 
functions in the two following subsections, reporting all the details of the computation in 
Appendix A and B. From these results we determine the FDR and the universal amplitude 
ratios TZr and TZc- 

A. The equation of motion and its solution 

At one-loop level the tadpole contribution T to the equation of motion (27) is 

m = J{dq)Cl{t,t) = 2^*dt'^|(dg)e^2^^(*-*') = 



fl-d/2 



6r(i 



(87r)ir(5 



= 2N,rA^-'^/^ 



(76) 



where Eq. (59) with r^ = has been used. For later convenience we introduce Nd 

2/[iAnr/'Tid/2)], 



'' = -Te^ 



Q S 

— --(7E + ln2) 
16 8^' ' 



+ 0{e) = ^ + co + 0{e), 
e 



and go = N^go- Therefore Eq. (27) becomes 



= dtin + -lit" + goTdt^ '^^'^m . 
o 



(77) 



(78) 



We look for a solution of the form m{t) = mG(t)[l + gop(t)] + 0((?q), where niQit) = 3/(2t) 
is the Gaussian solution. It is straightforward to show that p(t) = r^/iS — d/2)t'^^'^, so that 



™w = i/a 



1 - 9o: 



Td 



-f/2 



3-rf/2 
Expanding in e = 4 — li, we obtain the known result [14,15] 



0{9l) 



(79) 



m{t) 



1-5-0 

e 



l-5o (^Int + co-^) +0{e'^,ego,go' 



(80) 



As expected, a dimensional pole ~ 1/e appears in the expansion. To get finite renormalized 
quantities {M,g) in the limit e ^ one has to renormalize the bare quantities {Mo, go) 
according to the well-known renormalization procedure (see, e.g., Ref. [38]): Mq = Z^^'^M, 





1) 



(2) 



FIG. 3. One-loop diagrams contributing to the response function i?q(t,s). Directed lines rep- 
resent response propagators R'^{t,t'), whereas undirected ones represent correlations C^{t,t'). 
Causality implies that both diagrams vanish for s > t (as a result of the directed line joining 
the external legs of the diagrams). The corresponding expressions (reported in Appendix A) are 
indicated with Ii{t,s) and l2{t,s), respectively. 



go = ZgZ^^fi'^g where /i is an arbitrary momentum scale that we will neglect in the following, 
Z = l + 0{g'^), and Zg = l + 3g/(2e) +0{g'^). Taking into account Eq. (16) one recognizes in 
Eq. (80) the renormalization constant Zy^ jZ. Once the renormalized m(t) is expressed in 
terms of ^, its scaling behavior at the RG fixed point g* = 2e/3 + 0{e^) [for d > A, g* = and 
the predictions of the Gaussian model become exact] clearly shows up: m{t) ~ t~^ + 0{e^) 
with ^ = 1/2(1 + c_i^*) + 0(e2) = 1/2(1 - e/2) + 0{e^) = l3/{uz) [13] [we recall that 
z = 2 + 0(e2), l3/v=l- e/2 + ^2 = 1 - e/2 + 0{^)]. 



B. The response function 



At one-loop order the expression (58) for the response function (23) gets modified be- 
cause of the one-loop term contributing to the magnetization m{t), computed in Eq. (80). 
Accordingly, rri^if) = 3/(2t) - 3gordt^-l/{3 - f ) + 0{g^), and therefore 

<(t,s) = e-'^'(*-^)-/>*'™'(*') 



S\3/2 

~t) 



I +90 



t^ 2 ~ S 2 



3 - ^ 

"^ 2 



rd- 



+ 0{gi 



-c^it-s) 



(81) 



In addition to this contribution two further terms come from the interaction vertices of £i 
[see Eq. (19)] and are depicted in Fig. 3 [note that, up to this order, in computing the 
diagrammatic contributions one can use the tree-level expressions Eqs. (58) and (59) - with 
Tm = - for the response and correlation function]. In terms of these diagrams the response 
function reads 



i?q=o(t, s) = Rl^,{t, s) - |/i + goh + 0{gl) . 



(82) 



The relevant integrals have been calculated in Appendix A, Eqs. (A9) and (All). Collecting 
the three terms together we obtain, in terms of the renormalized coupling g {x = s/t) 



i?q=o(t,s)=a;3/2 1 



3, 3 /tt^ 7 ^ , , 

8 4 I 3 2 •'"^ ^ 



+ 0(^^e^,e 



(83) 
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where 

fni^) = 3 + ^ + 3^^^ ln(l - x) - 2Li2(x) , (84) 

2 X 

[we recall the definition of the dilogarithmic function Li2(a;) = Y^'^=ix"'/n'^] is a mono- 
tonically increasing function with /r(0) = and /ij(l) = 7/2 — 7r^/3. Note that all the 
dimensional poles and the non scaling terms originally present in the expressions of Ji^2 
cancel out in the sum, as they should. 

At the fixed point g* = 2e/3 + 0(e^) [38], the response function can be written as Eq. (34) 
with 

a = 0{e^) , (85) 

= -l + lr + Oie^) = -l + le + Oie^), (86) 

^« = l + r(^-|)+0(6^) = l + e(^-l)+0(e^), (87) 

fnix) = 1 + ^rfnix) + 0(e2) = 1 + i/^(x) + O(e^) . (88) 

6 agrees with the one-loop expression of —l35/{vz) [38]. Note that, as a difference with the 
case rriQ = [7], the response function present a one-loop correction to the Gaussian result. 
In Ref. [7] it was found that Ar = 1 + O(e^), leading to 

A first numerical estimate of TZr in three and two spatial dimensions can be obtained by 
setting e = 1 and e = 2 in the previous equation, giving TlR^d = 3) ~ 0.9 and IZR^d = 2) ~ 
0.8. However, sizable systematic shifts of these estimates might still come from higher-order 
terms in the e-expansion. 



C. The connected correlation function 

In terms of the Feynman diagrams depicted in Fig. 4, the connected correlation function 
at one-loop reads 

Cq=o(t, s) = Cl^oit, s) + |/5 + 9oh - f/r + 0{gl) . (90) 

Using the results reported in Appendix B and summing all the contributions we obtain 
the one-loop correlation function 



C^=o{t,s) = ^sx^/' 



^ ^ ^ ^ "320 ^ Teo'"' "81^^)+ ~9fci^) + C'if, e^, e'^ 



(91) 



where fc{x) is given by Eq. (B31) with /c(0) = 0. Eq. (91) can be cast in the scaling form 
(35), with e = -(36/{uz), a = 0{e^), 
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FIG. 4. One-loop diagrams contributing to the correlation function (see also the caption of 
Fig. 3). The diagrams represented here are invariant under the exchange s ^^ t, and therefore 
one can assume t > s. Their corresponding expressions are given, from top to bottom, by I^^t, s), 
lQ{t, s), and Ij{t, s) reported in Appendix B. 
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Taking into account that, up to one loop, Ac = 2(1 + e/6) + O(e^) [7], one finds 
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leading to the numerical estimates TZc{d = 3) ~ 0.21 and TZc{d = 2) ~ 0.18. 

D. One-loop FDR 

It is now easy to compute from Eqs. (87) and (92) the FDR 



A 



R 



TT 



2 \ 



Acil 



5 V600 looy ^ 



(92) 



(93) 



(94) 



(95) 



Let us comment on the main features of X°°. Close to rf = 4 it decreases as the dimensionality 
decreases. Probably a monotonically decreasing behavior has to be expected down to d = 1 
(the lower critical dimension) where it is known that X°° = 1/2 [37]. One can provide 
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numerical estimates of X°° in physical dimensions d = 3 and d = 2 hj evaluating Eq. (95) 
for e = 1 and e = 2, respectively, leading to X°°{d = 3) ~ 0.78 and X°^(d = 2) ~ 
0.75. These values are mainly indications of the actual ones, since higher-order terms in 
the e-expansion might give relevant contributions. However, the smallness of the one-loop 
correction (which is just the 3% and 6% in d = 3 and 2, respectively, of the mean-field result) 
strongly suggests that such estimates can work effectively, despite the low order considered 
in perturbation theory. For this reason in the next section we proceed to a numerical 
evaluation of X°^ by means of an extensive Monte Carlo simulation of the two-dimensional 
Ising model with Glauber dynamics. In Ref. [30] a phenomenological evolution equation 
for the fluctuating magnetization ((/^(q = 0,t)) [corresponding to a tree-level approximation 
of the actual Langevin equation (7)] is used [see Eq. (52) therein] to explain the values 
of X°° found within the Gaussian approximation. The resulting expression X°° = [2(3 + 
3vz)/{2f3 + Avz) for X°°(mo 7^ 0) can be naively used beyond the original approximation, 
giving X°° = {2(3 + 3i'z)/{2f3 + Auz) = 4/5 — e/50 + O(e^) which provides a rather good 
estimate of the actual result Eq. (95). 

VI. MONTE CARLO SIMULATIONS 

Monte Carlo (MC) simulations have been widely used to investigate the problem of the 
non-equilibrium critical dynamics and aging in two dimensions (see for instance [34,36,42-51] 
and references therein) providing results which are in rather good agreement with those 
obtained within different analytical approaches (see, e.g., Ref. [2]). It is then natural to test 
our predictions on the critical aging from a magnetized state with short-range correlations 
by comparing them with MC simulations of two-dimensional Ising model defined by the 
Hamiltonian (6) and with spin-flip dynamics. 

A. Details of the simulations 

We performed our Monte Carlo simulations using the Heat Bath updating rule, simulat- 
ing a large system oi N = L x L (with L = 10^ and periodic boundary conditions) spins Si, 
i = 1, . . . ,N and averaged our observables over more than 3000 different runs. We computed 
the magnetization 

M(t) = ^i^5.(t)^, (96) 

the connected two-time (auto) correlation function 

CA{ti,t2) = /^Y.S^itl)S,{t2)^ - M{ti)M{t2) , (97) 

where (. . .) stands for the MC average. We also considered the thermoremanent magneti- 
zation (TRM) at a site i and time t2 for a fleld applied at the same site from time t = to 
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p(ti,t2)= f'dt'RAit2,t'), (98) 

Jo 
where RA{t, s) is the (auto)response function of the system at time t and site i to a magnetic 
field hi apphed at the same site i but at an earher time s. In MC simulations a random 
probing field hi (with hi = and hihj = K^ 6ij) is usually used to extract the response 
function. However, measuring numerically the response in the case of a uniformly magnetized 
system is harder than in the usual case where M{t) = 0. Indeed one usually applies a small 
field hi at lattice site i, checks that the system remains in the linear regime, measures the 
local magnetization Mi = (Si) at a later time and finally deduces the response via Mi/ hi [49]. 
Here, however, one would need to compute [Mi{h) — Mi{0)]/hi which involves two different 
precise measures of the magnetization in the presence and in the absence of the field hi. 
Instead of using this numerically demanding procedure, we follow the approach proposed 
in Refs. [43,44], and later used in Ref. [47,48,52] which allows for the determination of the 
response function without introducing any magnetic field in the simulation. This approach 
has been adopted in Ref. [47] to compute X°° [see Eqs. (2) and (3)] for the two-dimensional 
Ising model with conserved and non-conserved dynamics leading to satisfactory results. For 
details on the method, we refer the reader to Refs. [44,47]. Here, we just recall the formula 
we used to compute the TRM: 

Tp{h,h) = j.Y.{SMAS,{U)), (99) 

i 

where the function ASi{ti) is computed during the simulation from time t = to ti and is 
defined by 

A5,(ti) = X: %), [S^is) - tanh(/ir(s)/T)] , (100) 

s=0 

where the function X(s) gives the index of the spin to be updated at time s (which depends on 
the random update performed in the MC simulation), T is the temperature (T = Tc in what 
follows) and hf{s) is the Weiss local magnetic field accounting for the effect of the interaction 
specified by the Hamiltonian (6) on the spin to be updated at time s (in the present case it 
is just the sum of the value of the neighboring spins, i.e., h^{s) = J2j:\i-j\=iSj{s)). Using 
this method, the functions 6*^(^1,^2) and p(ti,t2) can be calculated in the same simulation 
for any value of ti and t2, allowing for very efficient simulations. The initial state of each 
MC run is prepared such that {M(t = 0))ic = Mq (|Mo| < 1) where (. . .)ic stands for the 
average over the different realizations of the initial condition. In particular the value of each 
single spin in the lattice is chosen to be -|-1 and —1 with probability p^ = {1 + Mq)/2 and 
p_ = l—p+, respectively, leading to Aq = {[M(t = 0) — Mo]^)ic = {1 — Mq)/N, and vanishing 
initial correlations. Given that tq ~ Aq"^ [see Eq. (14)] one generally expects corrections 
to the leading scaling behavior which could in principle be reduced by preparing the initial 
state with sharply fixed magnetization. 

B. Scaling forms 

We now present the numerical results obtained using the method described in the pre- 
vious subsection. We first consider the scaling behavior of the magnetization and of the 
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FIG. 5. Scaling of the magnetization (left) and of the connected correlation function (right) for 
different initial conditions with magnetizations Mq. 

correlation as functions of the initial magnetization Mq. According to the theoretical pre- 
diction Eq. (30), we expect the magnetization to scale as 



r 



'FM{tM^ 



4.02 ^ 



101) 



where we have used the value of the critical exponent for the two-dimensional Ising model: 
(3 = 1/8, zy = 1, z = 2.1667(5), and 9 = 0.383(3) (see,e.g., Ref. [2]), leading to k ^ 0.249 
[see after Eq. (30)] and (3 /{vz) ^ 0.058. For the short-time behavior one expects Fm{x —>■ 



0) 



~ x" 



X 



0.249 



This scaling is rather well displayed by our data, as shown in Fig. 5. 



Having tested the initial increase of the magnetization, we now consider the connected 
correlation function. In the previous sections we focussed on the scaling properties of the 
correlation for vanishing momentum Cq=o(t, s), whereas now we are interested in the behav- 
ior of the (auto) correlation function CA{t,s) = /(dg)Cq(t, s). According to the discussion 
in Sec. Ill, the scaling form of Cq(t, s) is given by Eq. (33) where the scaling function has 
an additional dependence on the variable ?/ ~ g^(t — s). Accordingly we expect 



I3S 



CAit, s) = s{s/t)-{t - sr--Fcis/t, Mot^) 



(102) 



where the change a ^~>- a — d/z is due to the integration over the momenta. To make evident 
the scaling with respect to Mq, we consider the correlation Cyi(t,t/3) for different initial 
magnetizations Mq. Thus we expect a data collapse for t large enough according to 



-^^GcitM^) = t-°-iiGc(tMo'-°'). (103) 

■d/z = —2l3/{vz). The corresponding data are reported in 



CA(t,t/3) = t 

where we used the relation 1 + a - 

Fig. 5. Apart from a transient short-t effect, we observe, again, a rather good data collapse. 



C. Correlations, Responses and X°° 

We now specialize to completely ordered initial condition (Mq = 1). According to 
Eqs. (34) and (35), the two-time connected correlation function and the integrated response 
TRM are expected to share the same critical scaling form which we write as 
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FIG. 6. Scaling of the connected two-times correlation (left) and integrated response (right) 
functions for completely ordered initial conditions using Eq. (104). 



s'^/^^CA{t,s) = Ut/s) 



and 



.2^/^V(t, s) = Ut/s) 



(104) 



with 213/vz ^ 0.116 and fc,p{x > 1) ~ x'^ where (p = -a+f + ff = l+(3{6+2)/{uz) ~ 1.98 
(in two dimensions 6 = 15). This behavior is again rather clearly displayed by our data, as 
Fig. 6 shows. Note however that due to the fast decay of the correlations, it is difficult to 
obtain good data for x = t/s ^ 20 and that the TRM is affected by systematic deviations 
for the longer times. 

Finally, we consider the value of the FDR that can be alternatively read off from the 
limit 



X' 



hm^^ 

Ca^O CAit, s) 



(105) 



which indeed provides the estimate for X°° as defined in Eq. (3). The plot in Fig. 7 clearly 
shows that Tcp{t,s)/CA{t,s) in Fig. 7 is rather well approximated by a straight line in 
a wide range of x = s/t, both in the case Mq = and Mq ^ 0. (This point was al- 
ready observed in Ref. [36], see also [53]) The cross-over towards the quasi-equilibrium value 
Tcp{t, s)/CAit, s) = 1 takes place for x ~ 0.75-0.85 (although the crossover is not complete 
in the set of data for the case Mq = 0). A similar qualitative behavior is displayed by the 
field-theoretical Gaussian approximation for the same quantity, reported in Fig. 2. A more 
quantitative comparison, however, would require the analysis of the effects of fluctuations. 
In order to provide the numerical estimate for X°° we use only the set of data which gives 
reasonable errors bars (typically t/s ^ 20), reported in Fig. 7. From these data, we roughly 
estimate 






0.73(1) 



(106) 



which compares rather well with the one-loop prediction X°° ~ 0.75. Using the expression 
for X°° provided in Ref. [30] (see also Subsec. V D) as an improved Gaussian approximation 
one gets X°° = [3 + 2l3/{iyz)]/[A + 2(3/{vz)\ ~ 0.757 which is indeed rather close to the 
numerical result. 
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FIG. 7. Tcp{t, s)/CA_{t, s) versus s/t for disorder and non disorder initial conditions (left), 
our data we estimate X^q = 0.73(1) for ordered initial conditions (right). 



VII. CONCLUSIONS 

In this paper we have considered the non-equilibrium dynamics of a system in the Ising 
universahty class starting from a state with initial magnetization Mq and evolving at criti- 
cality according to a purely dissipative dynamics. By means of RG arguments we found that 



in the aging regime {t ^ s ^ t„ 



Mr 



-\JK 



) the zero-momentum response and connected 



correlation functions display the scaling behavior 

i?q=o(t, s) = AR{t- sf{t/sf^R{s/t) , 
Cq=o(t, s) = Ac s{t - s'nt/sffcis/t) 



(107) 
(108) 



independently of the actual value of Mq 7^ 0. We showed that the exponents 9 and 9 are 
not new non-equilibrium quantities, but they are given by standard equilibrium exponents 

as 

f35 



9 = 9 =-- 



vz 



(109) 



In addition we showed that the ratio of the non-universal constants Ar and Ac with the 
corresponding ones obtained in the case Mq = are universal. 

We solved exactly the dynamics within the Gaussian approximation and we derived 
all the relevant universal ratios and scaling functions. Then we considered the effect of 
fluctuations in order to provide more accurate predictions in physical dimensions c? = 2, 3. 
We computed the first order correction in e-expansion to the response and to the connected 
correlation functions and from them we obtained the universal ratios and associated scaling 
functions. In particular, for the FDR we found 



5 1 600 100/ ^ ^' 



;iio) 



that gives X°° ~ 0.78 in rf = 3 and X~ ~ 0.75 in rf = 2. On the other hand, if the FDR 
is defined from the non- connected correlation function then the asymptotic value would be 
zero because of the scaling in the volume and times of the additional term M(t)M(s). 
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In order to confirm our theoretical predictions we carried out extensive Monte Carlo 
simulation for the two-dimensional Ising model with Glauber dynamics using the approach 
of Refs. [43,44]. We checked the predicted scaling forms for the two-time response and 
correlation functions in addition to the well-known scaling behavior of the magnetization [16]. 
Our analysis leads to a quite accurate determination of the FDR 

X~c = 0.73(1), (HI) 

which is in rather good agreement with the field-theoretical estimate. 

It is worth noticing that at criticality X°° < 1 independently of the initial magnetization. 
Therefore, insisting on the definition of an effective temperature, one finds T^s = Tc/X'^ > 
Tc. This fact can be naively understood for the case of a quench from high temperature to 
the critical point: The system, because of its slow dynamics, is somehow unable to transfer 
properly heat to the thermal bath, resulting in an effective higher value of the temperature. 
In the case of a sudden heating from a low-temperature configuration (with non- vanishing Mq 
and given initial short-range correlations Aq) to the critical point one would have expected 
(as observed in some cases [52]) T^s being larger than the initial temperature but smaller 
than Tc. This is not the case here, showing that the properties of T^s at criticality are both 
non trivial and counter-intuitive. 

Let us now compare our results with those available in the literature. To our knowledge 
the ageing properties of the non-equilibrium critical relaxation from an ordered (alternatively 
magnetized) state has been the subject of three recent studies [31,30,9]. In Ref. [30] the 
lattice Ising model with Glauber dynamics has been solved in the mean-field cases of (a) 
fully connected lattice (infinite interaction range) and (b) nearest-neighbor interaction in 
the limit of large dimensionality d. The corresponding results are completely equivalent 
to our Gaussian solution, as expected. In Ref. [9] the ageing properties of the spherical 
model with relaxational dynamics have been considered. The scaling forms for response and 
correlation functions found there agree with our general predictions. In this case X°° is a 
monotonic function of the dimensionality d, interpolating between 1/2 at d = 2 and 4/5 at 
d = 4. This shares with our result for d < 4 a downward correction to the Gaussian result 
4/5. However a more quantitative comparison is clearly not possible, since the two models 
belong to different universality classes. The non-linear o"-model has been used to study the 
critical aging dynamics of ferromagnetic systems with 0{N) symmetry (A^ > 2) in Ref. [31]. 
Correlation and response functions are computed in a dimensional e-expansion around the 
lower critical dimension d = 2, with e = d — 2, leading to X°° = 1/2 + 0(e). These results 
do not compare directly with ours, since the behavior of 0{N) models is affected by the 
massless fiuctuation modes in the directions transverse to the decaying magnetization, in 
addition to the longitudinal one. In order to make a contact with Refs. [31,9] we have 
also studied the non-equilibrium aging dynamics of the 0{N) model up to first-order in the 
e-expansion. The analysis will appear elsewhere [54]. 
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APPENDIX A: FEYNMAN DIAGRAMS FOR THE RESPONSE FUNCTION 

There are two one-particle irreducible (IPI, i.e., that cannot be disconnected by cutting 
one propagator) diagrams entering in the calculation of the response function [see Fig. 3]. 
One is I{t) [see Eq. (76), depicted in Fig. 3(1)] and the other is the "bubble" in Fig. 3(2): 

BRcih, si) = Jidq)Rl{h, Si)C°(ti, si) . (Al) 

For ti > si [because of causality, BRc{ti < si, si) = 0] 

Bncih,s,)=l{dq) (^£i)'^%-'i^(*i-0^^e-'''(*^+^^) |J' dxxV^'^ (A2) 

= 2(87r)-iiJ,(ti,si) (A3) 

with 



HdUi, si) = tr"^ / dx x'^iti — x) 
Jo 



A,t^'^'/' - A,t^\h - s^y/' + B,r,\s,{h - s,r'^ + Cdtfsl{h - s,r'^ 



where 



+D,t^M(tl-Sl)-^+^/^ (A4) 

^^ " (2 + e/2)(l + e/2)e/2(-l + e/2) ' ^^^^ 

6 

^' ^ "(2 + e/2)(l + e/2)(-l + e/2) ' ^^^^ 

3 

^^ = ~(2 + e/2)(-l + e/2) ' ^^^^ 

Dd = --, rr- (A8) 

"^ (-l + e/2) ^ ' 

This parameterization is convenient because all the divergences arises only from those terms 
containing Ad or Dd- 

Thus the integrals for the response function read (as usual x = s/t) 



hit,s) = fdt'R%it,t')Iit')Rl^,it',s) = fdt'{t'/tf'\s/t'f''2NdTdt"-"' 

J s J s 

Inx+'eln^x + Oie^) , (A9) 

1/ ^ L 4 J 

corresponding to Fig. 3(1) and 



2iV,r,(s/t)3/2L__i__ = 2Ndrd{s/tf/h^/' 
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.(t,s)= / dti f ' dsiRl^Q{t,t,)V2mGiti)B{ti,si)V2mGisi)Rl^Qis,,s) 

J S J S 






(AlO) 



corresponding to Fig. 3(2). Using the parameterization (A4) for Hd-, I2 is reduced to five 
simple integrals whose sum is 






2, l-x vr^ /I \ 3, 
— ma; 1 h 3 ln(l — x) I f) H — In a; 

C Z: O \X / Z: 

+-ln2x-2Li2(x) +0(e) 



(All) 



with X = s/t and Li2(a;) = E^i^^^/n^. 



APPENDIX B: FEYNMAN DIAGRAMS FOR THE CORRELATION FUNCTION 



There is only one new IPI diagram for the correlation functions, depicted in Fig. 4(5) 
[we define t< = min{t,t'}]: 



^ ' ' dti dt2- 






{sny/^ Jo ^" Jo ^'^" {t + r-h- hy/^ ■ ^^^^ 

This enter in the connected diagram represented in Fig. 4(5) [we assume t > s, given that 

hit, s) = y dn J dT2 R%i^, n)V2mG{ri)Bcc{ri, T2)V2mG{T2)R%it^ ^2) 

dn r dr2{s/n)-'/'{3/ny/^BcG{ri,r2K3/r2f\t/r2)-'/' 
Jo 

+ Tdn f'dT2{s/T,)-'/\3/ny/'BcG{n,T2){3/T2Y/\t/T2)-'/' 
Jo Js 

= 3(st)-3/2 2 / dn r dT2TiT2BcG{ri,T2) + I dr^ I dT2TiT2Bcc{ri,T2) 
. Jo Jo Jo Js 

= 3{st)-^/^[2B + A] . (B2) 

It is easy to realize that both A and B do not diverge for d -^ 4. Therefore no dimensional 
regularization is required and one can calculate them directly in rf = 4: 
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where 
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4 Jo 



dyi {y^x 
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JO Jo 1 



3^3 



Jj 1 Jj 
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+ 1// -Xi -Xs)^ ' 



(B5) 



is a negative and monotonically decreasing function with A{x) = — x^/252 + 0(x'^) and 



^(1) 



1153 I 7r2 24 



ln2 ~ —0.016. [Although ^(x) can be expressed in terms of elementary 



3600 ' 30 25 

functions we do not report here its lengthy form.] Thus 



h{t,s) = 3{st)-'^^^[A + 2B] 



Sis 
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There are two diagrams involving Bjic{ti, ^i) [see Fig. 4(6)], namely 
h{t,s)=J dt^j^ dtiRl^Q{t,ti)V2mG{ti)Bjic{ti,t2)V2mG{t2)C^^=oit2,s) + 

y dti ^' dt2 Cl^oit, h)V2mG{ti)BRc{t2, h)V2mG{t2)R%{s, h) , (B7) 

that can be written as 



/6(t,s) = |(5t)-3/2 / ^^^^2 / dtitiBRc{ti,t2) + S^ dt2tf dtihBRc{tl,t 

UO Jt2 Js Jt2 

+ r dtiti r dt2t2BRcit2,ti] 

Jo Jti 



(B8) 



Using the parameterization of Eq. (A4), each of the three integrals decomposes in five 
integrals that can be easily computed. Therefore 



h{t,s) = l{str'/'[{A) + {B) + {C)] = ^-^^ E E( 
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where 
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Summing up one finds 
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where 



„, , 2x 7 , ,_, , oT. / X 31 41n l-x 1 
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is a monotonically increasing function witli B{x) = a;"^/63 + 0{x^) and B{1) = ffn ^ T 
0.068. 

Tlie diagram involving I{t) is [see Fig. 4(7)] 

hit, s) = 1^ dt' Rl^,it', s)Iit')C%it, t') + 1^ dt' i?°=o(t', t)/(t')CS=o(^', s) = 
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Finally there is the contribution coming from the one-loop correction to the magnetiza- 
tion. At this order in perturbation theory it can be simply evaluated as 



(B28) 



Using the value of m^(t) given by Eq. (79), we obtain 
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Summing up all the contributions according to Eq. (90) one finds 
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[/c(0) = 0] where A{x) and B{x) are given by Eqs. (B5) and (B26). 
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